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We examine fractional charge and spin errors in self-consistent Green’s function theory within a 
second-order approximation (GF2). For GF2 it is known that the summation of diagrams resulting 
from the self-consistent solution of the Dyson equation removes the divergences pathological to 
second-order Mpller-Plesset theory (MP2) for strong correlations. In the language often used in 
density functional theory contexts, this means GF2 has a greatly reduced fractional spin error 
relative to MP2. The natural question then is what effect, if any, does the Dyson summation have 
on the fractional charge error in GF2? To this end we generalize our previous implementation of 
GF2 to open-shell systems and analyze its fractional spin and charge errors. We find that like 
MP2, GF2 possesses only a very small fractional charge error, and consequently little many electron 
self-interaction error. This shows that GF2 improves on the critical failings of MP2, but without 
altering the positive features that make it desirable. Furthermore, we find that GF2 has both less 
fractional charge and fractional spin errors than typical hybrid density functionals as well as random 
phase approximation with exchange. 


I. INTRODUCTION 

Self-consistent single-particle electronic structure 
methods are of great interest because they combine 
conceptual and computational simplicity while being 
free of a starting reference bias. Probably the most 
well known example of this is density functional theory 
(DFT)m. Often these methods are designed to satisfy 
known constraints that an exact electronic structure 
theory should obey. For example, the total electronic 
energy should vary linearly in the fractional electron 
occupancy between integer electron numbers [2t{5], while 
for a simple one-electron system like the hydrogen 
atom the energy should be degenerate with respect to 
variations in the fractional spin|6]. To the extent that 
a method disobeys these two exact constraints such a 
method will display many electron self-interaction and 
static correlation errors, respectively. To understand 
how these errors are connected to fractional electron 
behavior, let us consider the H 2 molecule in the infinite 
dissociation limit in two different scenarios: 

First let us focus on a fractional charge error. If an 
extra electron were placed on this system with fractional 
occupancies on each H atom given by n a , n&, n a +n 5 = 1, 
then the net energy change of the system would be 
Fa = — ( n a A + n^A) = — A , where A is the electron 
affinity of the hydrogen atom. The energy of each sub¬ 
system would vary linearly in the occupation, and the to¬ 
tal energy of the entire [H 2 ] - system would be the same 
regardless if the extra electron were delocalized across 
both atoms, or localized only on one. If the energy of the 
subsystem varied nonlinearly, then either the delocalized 
or localized solution would become unphysically lower in 
energy depending on whether the curve was convex or 
concave. This unphysical behavior would be a simple 
manifestation of many electron self-interaction error. 

Let us examine now a fractional spin error. For neu¬ 
tral singlet H 2 in the infinite dissociation limit each H 
atom should have half a spin up and down electron, and 
the energy of the singlet should be identical to that of 


broken-symmetry solutions where spin up and spin down 
electrons have localized on different atoms. Therefore 
a method’s failure to yield equivalent energies for one- 
electron H with fractional or integer spin is equivalent to 
the failure to describe multireference static correlation 
energy. 

From this simple example it is clear that fractional 
charge and fractional spin errors are deeply connected 
to many electron self-interaction and static correlation 
errors. For this reason these errors have been stud¬ 
ied extensively, and are known to have severe nega¬ 
tive consequences for a given method’s description of 
properties that depend on electron delocalization and 
static correlation effects [6LFTQ] . This language of frac¬ 
tional charge and spin has traditionally been used ex¬ 
clusively within the DFT community to analyze ap¬ 
proximate density functionals. However these concepts 
have begun to make inroads into other areas such as 
wavefnnction [IT] and many-body theory [T21IT6] . and den¬ 
sity matrix theory as wellflTj. For example it has been 
shown that MP2 (second order Mpller-PlessetjJE]) pos¬ 
sesses relatively little fractional charge error, but dis¬ 
plays a massive diverging fractional spin error (which 
is an alternative way of stating that MP2 diverges for 
strong correlations) [12] . Double-hybrid density func¬ 
tionals that include some MP2 correlation are expected 
to have a similar fractional electron behavior p33 [20] . 
The Random Phase Approximation fRPA] [2Tj on the 
other hand has minimal fractional spin error, but a se¬ 
vere fractional charge error [T3]. The closely related GW 
approximation [22] displays a fractional charge error sim¬ 
ilar to RPA, but a larger fractional spin error [16]. 

Motivated by these works, we find it interesting to 
extend this analysis to self-consistent Green’s function 
theory in a second order approximation (GF2)[23, [24]. 
Similar to MP2, GF2 includes all diagrams to second or¬ 
der in the bare electron-electron interaction, as shown in 
Figure [l] and is therefore exact for one electron systems. 
However in contrast to MP2 these diagrams are evalu¬ 
ated with self-consistent Green’s functions, obtained by 
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FIG. 1. Feynman diagrams for the second order self-energy in 
GF2. Here a red wavy line represents a two electron integral, 
while a black arrow line represents a Green’s function. From 
left to right the diagrams shown are the first order Hartree 
and exchange diagrams, and the second order pair bubble and 
second order exchange. 


greater fractional spin error. Any hybrid of these two 
fB3LYPf28l[3Q. or PBEh |29ll3T] . for example) will essen¬ 
tially trade one error for the other to the extent that the 
HF-type exchange is admixed in place of DFT exchange. 
What is worth noting is that, in the language of hybrid 
DFT, the Fock matrix in GF2 contains full HF-type ex¬ 
change (which in Green’s function theory is usually re¬ 
ferred to as first order exchange) yet we will show GF2 
yields both less fractional charge and fractional spin error 
than HF, B3LYP, and PBEh. This unique result comes 
about from a combination of the Dyson summation with 
including all diagrams to second order. 


iterative solution of the Dyson equation 


II. SPIN UNRESTRICTED GF2 THEORY 


G(uj) — G o(uj) + Go(cj)F(cj)Go(cj) 
+Go(<^’)F(cj)Go(cj)F(u;)Go(^’) + • • • 

= Gq(uj) (^(S(w)Go(w)) ) 

' n ' 

= [GoH " 1 - S(w )] -1 


( 1 ) 


To study open-shell systems we generalize G(u) to 
have two spin blocks 


G = 



0 

Gfl. 


( 2 ) 


Here Go(uj) is the Green’s function of a non-interacting 
system, while E(u;) is the proper self-energy, which in 
GF2 is truncated at second order and written as an ap¬ 
proximate functional of the Green’s function, E[G(u;)]. 
Because of the structure of the Dyson equation, the self- 
consistent G(uS) will contain an infinite order summation 
of the second order proper self-energy parts, E(cj). As we 
recently showed, this summation of diagrams allows GF2 
to give reasonably fine results for strongly correlated sys¬ 
tems such as stretched hydrogen lattices [23 when MP2 
would diverge. In the language of fractional electron er¬ 
rors, this suggests that GF2 improves tremendously over 
MP2 for fractional spins as a result of the self-consistent 
infinite order summation. An interesting question that 
arises then is what effect does this Dyson summation 
have on the more general fractional electron behavior? 
Relative to other methods such as RPA, GW, approxi¬ 
mate DFT, and Hartree-Fock (HF), MP2 has only a very 
small fractional charge error [12], and consequently lit¬ 
tle many electron self-interaction error (SIE). Ideally one 
would hope that GF2 improves on the disastrous frac¬ 
tional spin error of MP2 without deteriorating MP2’s 
impressively small fractional charge error. To investi¬ 
gate this question, here we will generalize our previous 
GF2 implementation [23] to open-shell systems and then 
investigate its fractional charge and spin behavior. 

Before closing this section it should be emphasized that 
what is challenging about the fractional charge and frac¬ 
tional spin errors is that any attempt to reduce one er¬ 
ror tends to exacerbate the other [13] l25l 126] , For exam¬ 
ple, a semilocal DFT functional (such as BLYP [27l [28] . 
or PBE|29j) will tend to have a large fractional charge 
error but a relatively smaller fractional spin error. On 
the other end of the extreme Hartree-Fock will have 
significantly less fractional charge error but a much 


where the spin-up and spin-down blocks are given by 

G„(w) =[(m <t +w)S-F <j -S <t H]“ 1 , o = a,p (3) 

The off-diagonal spin-blocks of G(cj) here are identically 
0 , meaning we do not allow for the possibility of spin- 
flips, and our solutions are constrained to be eigenstates 
of S z . In Eq. [ 3 ] S and F a are the overlap and Fock ma¬ 
trices, E cr (cj) is the self-energy, fi a is the chemical po¬ 
tential, and uo is an imaginary frequency. By introduc¬ 
ing /i a , /ip as separate chemical potentials we can al¬ 
low for different numbers of electrons in the respective 
correlated density matrices, P a , P^, which are given 
by Per = — G a (r=i/k B T), a = a,/3, where G cr (r) is 
the Green’s function fast Fourier transformed (FFT) to 
the imaginary time domain, and 1/^T is the inverse- 
temperature. The expression for ¥ a is the standard re¬ 
sult from spin-unrestricted HF theory, 

Ftj = hij + Y,( p ki + - PkiVikij , 

H (4) 

K = % + + p k>m - Pfrmj . 

kl 

However, unlike HF theory the density-matrices that en¬ 
ter this expression are those obtained from the Green’s 
function and thus include electron correlation effects from 
solving the Dyson equation. This covers the electron- 
electron interaction from zeroth through first order (the 
first order diagrams in Figure [l] are described by the HF 
mean-field). At second order in GF2 the electron-electron 
interaction is described by the frequency dependent self¬ 
energy, which is given in the imaginary time domain as 
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S “j( T )= X - G mn( T ) G kl( T ) G pq(- T ) v imqk{vipnj - V np lj) 
klmnpq 

- G mn( T ) G ll( T ) G p q (- T )VimqkVlpnj , 

^ij( T ) = y ' — G mn(. T ) G kl( T ) G pq( — T ') Y i r nqk{ v lpnj — Vnplj) 
klmnpq 

- G mn( T ) G kl( T ) G %q(- T )Vim q kVlpnj , 


(5) 


The reasoning for this expression is that at second or¬ 
der a spin-up (down) electron can have a pair-bubble 
interaction with both spin-up and down electrons, yet 
the second-order exchange can only proceed between like- 
spin electrons because we do not allow for the possibil¬ 
ity of spin-flips. Once X a (r) has been built it can be 
FFT to the frequency domain, and then we can rebuild 
Gcrfo;) with Eq. [3j Looking at Equations |3j [4j and |5j 
the spin-up and spin-down Green’s functions are coupled 
by the fact that Gp appears in the expression for 
(likewise with Gq, for X^), and that Gp contributes to 
F a through P/3 = —Gp{r=i/k B T) (likewise with G a , F a 
and F/ 3 ). Therefore Equations [3j [4] and [ 5 ] will need to be 
solved self-consistently, and at every iteration fi a and /i /3 
will need to be adjusted to give the desired number of a 
and (3 electrons. To start this self-consistent procedure 
we use a HF Green’s function generated by output from 
the Dalton program [32]. Note when G a = G/3, then 
= SI /3 and Eq. [5] reduces to the familiar expression 
for spin-restricted GF2[23|24]. 

The energy is evaluated as 


E = -Tr [(h + F a )F a + (h + F^P/ 3 ] 


Re Tr[G a (cj n )X) Q; (cj n ) + Gp(u n )^p(w n )\ 


(6) 


where uj n is a Matsubara frequency, uj n = ( 2 n + 1 ) 7 t&bT. 
This can be understood as essentially a spin-unrestricted 
HF-like energy expression supplemented with a frequency 
dependent correlation contribution from GF 2 . However 
it should be understood that all quantities are evaluated 
using the correlated P a as obtained from G a . 

Because the GF 2 approximation includes in the proper 
self-energy all exchange and Coulomb type diagrams to 
second order, it is by construction exact for one-electron 
systems, i.e. it is one-electron self-interaction and self¬ 
correlation free. Less clear however is its many electron 
self-interaction error (SIE) for general systems. To this 
end in the following we investigate the fractional electron 
behavior of GF 2 for several archetypical cases, and com¬ 
pare to standard density functional theory calculations 
ran with Gaussian 09[33j. 


III. RESULTS 


First we consider a single hydrogen atom with frac¬ 
tional spin up and down electron occupations, n a , up, 
that are varied in the interval 0 < n a < 1.0 and 
0 < ri /3 < 1.0. For this case with an exact method the 
energy should change linearly in the fractional electron 
number n = n a + np, while for constant n it should be 
invariant with respect to changes in the fractional spin 
m = n a~ n i3- Furthermore, there should be a discontinu¬ 
ity in the slope dE(n)/dn across the line n a + np = 1.0. 
Hence the resulting energy surface should be two flat 
planes that intersect along a seam. In Figure [2] we show 
the exact result, compared against that obtained with 
spin-unrestricted GF2. If first we focus only on the edge 
of the plane where one occupation is held fixed at inte¬ 
ger values of 0 or 1.0, this corresponds to the fractional 
charge behavior of a method and thus SIE. It is clear that 
GF2 reproduces the exact linear behavior almost per¬ 
fectly, similar to the result for MP2 in Ref. [12]. Now if 
we focus on the more interesting region towards the inte¬ 
rior, discrepancies between GF2 and the exact flat-plane 
behavior become apparent. To see this more clearly, in 
Figure[3]we plot the difference between GF2 and its inter¬ 
polated flat-plane surface. Looking at the diagonal region 
connecting the coordinates {n a = 1.0, np = 0.0} and 
{n a = 0.0, Tip = 1.0} we find a hill of fractional spin er¬ 
ror, where GF2 is not able to fully recover the static cor¬ 
relation energy. On either side of the hill ( n a +np < 1.0 
and tIq, T Tip > 1.0) we find shallow valleys where GF2 
moderately overestimates the correlation energy. The 
GF2 results here for fractional spin are in severe con¬ 
trast to the MP2 result from Figure 4 of Ref. [Tl]: which 
rapidly diverges to — oc correlation energy as one moves 
towards the center at {n a = 0.5, np = 0.5}. 

For purpose of comparison, we also show in Figure [3] 
the surface obtained with PBEh, as well as the GF2 re¬ 
sult obtained when the second-order exchange (SOX) di¬ 
agram is neglected in the self-energy (we call this GF2- 
NoSOX for short). Neglect of the SOX diagram intro¬ 
duces an unphysical one electron self-correlation error 
at second order. Because many electron SIE and static 
correlation are known to be connected^3] [26] [341436] . 
GF2-NoSOX should provide an interesting contrast to 
standard GF2. First let us compare the energy land¬ 
scapes for GF2 and PBEh. Despite being a hybrid func- 
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FIG. 2. The total energy for hydrogen with fractional spin up 
and down occupations n a , np evaluated with aug-cc-pVDZ. 
Top: the exact result. Bottom: GF2. 


tional, PBEh still yields a significant fractional charge 
error along the outer rim of the surface, along with a 
pronounced hill of fractional spin error running through 
the center. The GF2 energy landscape in contrast is ap¬ 
preciably flatter, with less fractional spin and fractional 
charge error. This is notable, because it has been stressed 
that simultaneously reducing the fractional charge and 
spin errors in a single-particle method is very difficult T31. 
For example, RPA+X (RPA with HF-type exchange) 
greatly improves over the fractional charge error of RPA, 
but at the price of gaining a considerably larger fractional 
spin error that is comparable to Hartrcc-Fock| 1 3 . Now 
let us examine the GF2-NoSOX result. As mentioned, 
neglect of the SOX diagram introduces a self-correlation 
error at second order, and as a result GF2-NoSOX gives 
significant fractional charge errors comparable to PBEh. 
At the same time, the hill of fractional spin error is ap¬ 
preciably reduced relative to both PBEh and GF2. 

To see the static correlation error more clearly, in Fig¬ 
ure [4] we plot the energy with respect to fractional spin 
m = n a — np for hydrogen with GF2, GF2-NoSOX, 
PBEh, B3LYP, and HF. This corresponds to the line 
running along the energy surface from the coordinates 
{n a = 1.0, np — 0.0} to {n a = 0.0, np = 1.0}. As 
stated previously an exact method should give a flat en¬ 
ergy curve from m = — 1 to m = 1. For example, HF 
has a massive hill at m = 0, which reflects the com¬ 
plete absence of static correlation energy in this method. 
Hybrid DFT yields much smaller fractional spin errors, 
with B3LYP being slightly lower than PBEh, likely be¬ 
cause it includes less HF-type exchange than PBEh. In 


comparison the fractional spin error is relatively lower 
for GF2, as a result of the infinite order summation from 
the Dyson equation recovering some static correlation. 
GF2-NoSOX has a much reduced fractional spin error 
relative to all four methods, with its energy at m = 0 
being not much different from that at m = ±1. It has 
been understood for some time that SIE can mimic static 
correlation p~3l [26[ l34lf36] . Usually this is considered in 
the context of SIE resulting from incomplete cancella¬ 
tion of Coulomb and exchange terms at first order. GF2- 
NoSOX’s small fractional spin error in contrast is purely 
arising from incomplete cancellation of Coulomb and ex¬ 
change terms at second order, resulting in an unphysical 
one-electron self-correlation that mimics static correla¬ 
tion energy, analogous to the situation that occurs with 
RPA+SOSEX|26] (second order screened exchange). We 
think this is an example of getting the right result for the 
wrong reason. Furthermore, GF2-NoSOX must obtain 
this slightly reduced fractional spin error at the price of 
gaining a tremendous fractional charge error, which is not 
a desirable trade. In contrast, from comparing Figures [4] 
and [3] it is clear that GF2 has less static-correlation error 
than typical hybrid density functionals, and importantly 
achieves this while being essentially one and two electron 
self-interaction free. This means GF2 genuinely recovers 
some static correlation energy, rather than fortuitously 
exploiting spurious self-interaction or self-correlation. It 
is worth mentioning that a similar analysis has been per¬ 
formed for RPA and RPA+X[f3l. Comparing our GF2 
result in Figure [4] to Figure 1 in Ref. 113], GF2 yields an 
appreciably smaller fractional spin error than RPA+X, 
but is still larger than RPA. 

Finally, to analyze the fractional charge error in closer 
detail we consider the energy of helium with respect 
to fractional electron number. In Figure [5] we plot 
E x (n ) — E x n , where E x (n) is the energy from method 
X for electron number 1.0 < n < 2.0, while E x n is the 
linear-interpolation between integer points with the same 
method. A method is said to be M-electron SIE free if 
E x (n) — E x n = 0forM — l<n< M\ 4]. We find that 
GF2 has a very small concave curvature. In contrast HF 
is moderately concave, while PBEh, B3LYP, and GF2- 
NoSOX are significantly convex. This clearly establishes 
that GF2 is one and almost perfectly two electron SIE 
free. Interestingly, comparing Figure [5] to Figure 4 in 
Ref. 43]- it is reasonable to conclude that GF2 should 
have less SIE than both RPA and RPA+X. 


IV. CONCLUSIONS 

We have analyzed fractional electron errors in self- 
consistent Green’s function theory by generalizing our 
previous GF2 implementation [23] to open-shell systems. 
Overall we find that GF2 has a very small fractional 
charge error, and a moderate fractional spin error. In 
comparison to other well known methods, we find that 
GF2 has both less static correlation and self-interaction 
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FIG. 3. The energy difference A E = E — Ey in for the hydro¬ 
gen atom with fractional spin up and down occupations n a , 
ri/ 3 , where E is the energy evaluated with fractional electron 
number, and Eu n is the flat-plane linear interpolation for a) 
GF2, b) GF2-NoSOX, and c) hybrid PBEh. All calculations 
are with aug-cc-pVDZ 


error than hybrid density functionals such as B3LYP and 
PBEh, as well as RPA+X and HF. Because the GW ap¬ 
proximation is diagranimatically identical to RPA, GF2 
will very likely have significantly less self-interaction er¬ 
ror than GW as well. Furthermore, it has been shown 
that CCSD has a roughly similar fractional charge error 
to MP2 m- From this it stands to reason that GF2 and 
CCSD will have comparable fractional charge errors. 

Essentially, by virtue of the Dyson summation GF2 
greatly improves on the tremendous fractional spin er¬ 
ror of MP2, but without deteriorating MP2’s relatively 
excellent fractional charge behavior. These results could 
suggest a way towards removing the fractional spin error 
from double hybrid density functionals |19J [20'. As a fur¬ 
ther salient point, GF2 is fully self-consistent and thus 
the converged density should reflect the relative lack of 
many-electron self-interaction error in a second order ap¬ 
proximation. MP2 in contrast is by definition a pertur¬ 
bative scheme that does not revise the underlying mean- 
field reference, and thus inherits the Hartree-Fock den- 



FIG. 4. The energy difference A E = E x — E x t for the frac¬ 
tional spin hydrogen atom, where E x is the energy from 
method X evaluated with the fractional spin population 
m = n a — np, n a + np — 1, and E x t is the energy from X 
with integer m — El. All calculations are with aug-cc-pVDZ. 



FIG. 5. The energy difference A E = E x —E x n for the helium 
atom, where E x is the energy from method X evaluated with 
fractional electron number, and E x n is the linear interpolation 
from method X. All calculations are with cc-pVTZ. 
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sity with its bias towards localization. This suggests that 
GF2 could find good application for properties that sen¬ 
sitively depend on electron delocalization. For example, 
much of the interesting physics in transition metal com¬ 
plexes depends on the slight delocalization of unpaired 
d electrons onto ligands, which is determined by the in¬ 
terplay of dynamic correlation and self-interaction error 
effects [37], 

Conceptually, GF2 is essentially a self-consistent 
single-particle theory where the energy is expressed as 
a functional of the single-particle Green’s function E[G\, 
in obvious analogy to DFT with density functionals E[p\. 


In terms of the energy, it is fair to say GF2 is a “Green’s 
function functional” with desirable fundamental proper¬ 
ties compared to standard hybrid density functionals. 
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